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Abstract 

The dielectric properties of single-shell spherical cells with an intrinsic dielec- 
tric dispersion has been investigated. By means of the dielectric dispersion 
spectral representation (DDSR) for the Clausius-Mossotti (CM) factor, we ex- 
press the dispersion strengths as well as the characteristic frequencies of the 
CM factor analytically in terms of the parameters of the cell model. These 
analytic expressions enable us to assess the influence of various model param- 
eters on the electrokinetics of cells. Various interesting behaviours have been 
reported. We extend our considerations to a more realistic cell model with 
a graded core, which can have spatial gradients in the conductivity and/or 
permittivity. To this end, we address the effects of a graded profile in a 
small-gradient expansion in the framework of DDSR. 

PACS Number(s): 87.18.-h, 82.70.-y, 77.22. Gm, 77.84.Nh 
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I. INTRODUCTION 



The interaction of the polarization of biological cells with the applied fields has resulted in 
a wide range of practical applications from manipulation, trapping to separation of biological 
cells [1], and even nanotechnology [2]. When a biological cell in medium is exposed to an 
applied electric field, there is an accumulation of charge at the interfaces and hence a dipole 
moment is induced in the cell. The strength of the polarization depends on the frequency 
of the applied field as well as on the permittivities and conductivities of cells and medium. 
The situation becomes more complicated when we consider structured particles because 
biological cells are usually modeled as conductive spheres (cytosol) with a thin insulating 
outer shell (membrane), assuming the shell is an isotropic, non-dispersive dielectric with 
conductive losses. In this case, there are additional frequency-dependent changes in the 
polarization. 

The Clausius-Mossotti (CM) factor determines the polarization of a biological particle 
in a surrounding medium, and is a measure of the dielectric contrast between the parti- 
cle and medium. The CM factor is important in biophysical research because it is closely 
related to the alternating current (ac) electrokinetic behaviors of biological cells, namely, 
dielectrophoresis [3], electrorotation [4], electro-orientation [5], electrofusion [6], as well as 
electrodeformation [7]. Any change in the cell's properties such as the mobile charges, or 
particle shape as well as the variation of medium conductivity or medium permittivity will 
change the CM factor, which is in turn reflected in the ac electrokinetic spectra. These spec- 
tra show characteristic frequency-dependent changes amongst other complicated features. 

Moreover, the conductivities and permittivities can have characteristic frequency depen- 
dencies due to the presence of mobile charges in membrane. Thus the constancy of these 
quantities is only an approximation and these quantities do change with frequency, giv- 
ing rise to additional dispersions. In this work, we aim to establish a dielectric dispersion 
spectral representation (DDSR) for the single-shell spherical cell model with an intrinsic 
dielectric dispersion in the cytosol. The DDSR was pioneered by Maxwell [8] in 1891 in the 
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context of interfacial polarization. When two media are put in contact (thus forming an 
interface) and an electric field is applied, polarization charge is induced at the interface due 
to the dielectric contrast between the two media. Although Maxwell considered a two phase 
system in which one phase is insulating, it can be readily generalized to a more general case 
when both media have complex dielectric permittivities. 

The DDSR was subsequently extended to spherical particles by Lei et al. [9] and further 
elaborated by Gao et al. [10] for cell models without shells, the single-shell model has been 
widely used to mimic a living biological cell as a homogeneous, nondispersive spherical par- 
ticle surrounded by a thin shell corresponding to the plasma membrane. The DDSR enables 
us to express the CM factor analytically in terms of a series of sub-dispersions, each of which 
with analytic expressions for the dispersion strengths and their corresponding characteristic 
frequencies expressed in terms of the various parameters of the cell model [9,10]. Thus this 
representation enables us to assess in detail the influence of the various model parameters, 
including structural, material, as well as dynamic properties of cells, without the need to 
analyze the full dielectric dispersion spectrum. 

The paper is organized as follows. In the next section, we review the dielectric dispersion 
spectral representation (DDSR) for the CM factor of an unshelled spherical cell model [9] 
to establish notations. We express the dispersion strength as well as the characteristic 
frequency of the CM factor analytically in terms of the parameters of the cell model. Then 
an intrinsic dielectric dispersion is included in the cell [10]. In Section III, we analyze the 
single-shell model with a dispersive core and a non-dispersive, insulating shell. We apply 
DDSR to the CM factor to obtain the analytic expressions for the dispersion strengths and 
characteristic frequencies. These expressions enable us to assess the influence of various 
model parameters on the electrokinetics of cells. In Section IV, we examine the influence of 
the individual parameters, such as the conductivities of the external medium and the cytosol 
on the dispersion spectra. Various interesting behaviours will be obtained. In Section V, 
we extend our considerations to a graded core, namely, the core can have spatial gradients 
in the conductivity and/or permittivity. We address the effects of a graded profile in a 
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small-gradient expansion in the general framework of DDSR. Discussion and conclusion will 
be given in Section VI. 



II. DIELECTRIC DISPERSION SPECTRAL REPRESENTATION 

In this section, we review the dielectric dispersion spectral representation for the CM 
(Clausius-Mossotti) factor of an unshelled spherical cell model [9]. The dipole moment p of 
a single sphere in uniform electric field [11] 

e e U D 3 

V = E (1) 

where e e is the permittivity of the external medium, D is the diameter of the particle and 
E is the electric field strength. U is the CM factor due to the dielectric discontinuity and 
follows the equation 

U = (2) 
€i + 2e e 

where q is the permittivity of the particle. In AC applied fields, we replace the permittivities 
with their complex counterparts: 

-> e* = e { + — , (3) 
too 

e e ^e* e = e e + — , (4) 
ten 

where i = o"« and a e are conductivities of the particle and of the external medium 

respectively. Then 

U ^u* = 6 * 1 - e * e (5) 
e* + 2e* e 

This gives the dielectric relaxation of a single spherical particle 

U* = U + - Ae , (6) 

1 + iuj I UJ C 

with the characteristic frequency u c and dispersion strength Ae: 
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o-i + 2a e 

u c = ^ — , (7) 

Ae = ^"^ - ii^L. (8) 
(jj + 2(T e 6j + 2e e 

It is related to the Maxwell- Wagner structure relaxation u c = 10 4 s _1 ■ • • 10 9 s _1 . 
The angular velocity f2 of electrorotation is 

= _£fM/ m [/* ( 9 ) 

where rj is the coefficient of viscosity. Note that Im U* < gives co-field rotation while 
Im U* > gives anti-field rotation. 

Then, when an intrinsic dielectric dispersion is included in the cell [10], we again replace 
the permittivities with the appropriate complex counterparts: 

e i =e i + i T • / + — ' 10 ) 
1 + tuj/uj c tUJ 

< = e e + % (11) 

The corresponding complex CM factor U int * can then be expressed in the dispersion 
terms as 

£Ant* = C/int + E , (12) 

where Ui nt = (e* — e e )/ (e* + 2e e ), Ae„s are the dispersion strengths and u n s are the charac- 
teristic frequencies. 

To solve for the dispersion strengths and the characteristic frequencies, assume the sum- 
mation term in Eq. (12) is of the form 

TJ. * TJ. - P ° ± PlW ( U ) 

Umt ~ Umt ~ l + R lW + R 2 w* (13) 



Po + Piw 



(14) 



(1 + w/ui)(l + w/uo 2 ) 
where w = iui and, Ps and Rs are constants in terms of the model parameters. 
For the characteristic frequencies, solve the following quadratic equation 

1 + R x w + R 2 w 2 = (15) 



and the u n s are minus the solutions to the equation. They come out to be, in terms of the 
model parameters, 

001 = 9(9 ~ J 2 ^ + <T i+ ( Ae * + 2e « + + ^ ( 16 ) 

c^ 2 = , 1 J 2tr e + a, + (Ae, + 2e e + ei )w c - vf], (17) 
ZyZ€ e + 6j J 

where 

T = -4(2e e + e;)(2cr e + ai)uo c + [2cr e + ^ + (Ae 4 + 2e e + e^] 2 . (18) 

For the dispersion strengths, performing partial fraction can express the summation term 
in the form of the summation term in Eq. (12). The dispersion strengths turn out to be, in 
terms of model parameters and characteristic frequencies, 

_ 3(-6^o- e ^i + e e ajUJi + €ja e u c - e e (TiUJ c + AQe e aw c ) 

1_ (26 e + 6 l )V(^ 1 -C 2 ) ' 

^ e _ 3(e i (T e cc'2 — e e <7iUJ2 — ei<7 e uj c + e e <JiUj c — Aeje e u; 2 w c ) 
(2e e + ei) 2 u 2 (uji - u 2 ) 

It is worth remarking that, two dispersion terms appear in Eq. (12): the first term (i.e. 
when n — 1) is due to the phase difference between the cell and the medium, and the second 
term (i.e. when n = 2) is due to the presence of the intrinsic dispersion inside the cell. 

This is a special case of the model mentioned in the following section. It is interesting 
to compare that this model, with no shell, but the same core as the next model, has two 
dispersion strengths (and the same number of characteristic frequencies), while the next 
model, with shell, has three dispersion strengths (and the same number of characteristic 
frequencies) . 

Similar work was done by Foster et at. [12]. For the case of nondispersive particle and 
medium, our solutions are indeed equivalent to those of Foster et al. However, for the case of 
dispersive particle and nondispersive medium, we quoted the exact analytic solutions while 
Foster et al. only presented the approximate solutions obtained by expanding the exact 
solutions using Taylor's expansion (cf. Section b of Ref. [12]). 
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III. SINGLE-SHELL SPHERICAL CELL MODEL, WITH A DISPERSIVE CORE 



The CM (Clausius-Mossotti) factor of an isotropic model with a non-dispersive homoge- 
neous core has been investigated [13,14]. Here we would like to establish the DDSR (dielectric 
dispersion spectral representation) of an isotropic model with a dispersive homogeneous core 
covered with a non-dispersive, insulating membrane [15]. 

The idea of DDSR is to mathematically extract the analytic expressions of the dispersion 
strengths and the corresponding characteristic frequencies from the CM factor. The CM 
factor for a single-shell spherical cell with isotropic, lossless dielectric membrane is [13,14] 

u = (2e m + e,j)(e m - e e )R 3 e + - e m )(2e m + e e )Rf 
lso ~ (2e m + ei )(2e e + e m )Rl + 2(e, - e m ){e m - e e )Rf 1 j 

where e is permittivity and R the radius; the subscripts e, m and i correspond to the external 
medium, the membrane and the cytosol respectively. 

For adaptation to our concerned model, the real constants e e , e m and 6j are replaced by 
the complex counterparts 

< = +~ (22) 
4 = e m + — (23) 

< = e e + % (24) 

ioj 

The complex e* contains the dispersive term ( ) to acc °unt for the intrinsic dispersive 
nature of the cytosol, while both the membrane and the external medium are non-dispersive. 
The CM factor becomes complex and can be written as 



^ 1 + liOjLOt 



where Ae t is the dispersion strengths, and cu t is the characteristic frequencies. 

Ae t and uj t can be solved easily using Mathematica. Assume the summation part to be 
of the form 
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u* -u- - Bo + BlW + B2Vj2 m) 



1 + Aiiy + A 2 w 2 + A 3 w 3 
B + B lW + 5 2 w 2 



(27) 



(1 + w/ui)(l + w/uj 2 ){1 + w/lj 3 ) ' 

where w = iuj and the As and _Bs are constants in terms of the parameters of the model. 
Performing partial fraction can express this term in the form of the summation in Eq. (25). 
To solve for u t , solve the cubic equation 

1 + A t w + A 2 w 2 + Aw 3 = (28) 

u t s are minus the solutions to this equation. 
Aei in terms of the constants Bs and u t is 

= (Bp + ui(-Bi + B 2 ^i))^2^3 ^ 

The rest of the Ae t s follow by cyclic permutation of the variables, namely 1^2,2^ 
3 and 3^1. 



IV. THE INFLUENCE OF INDIVIDUAL PARAMETERS 

This model depends on the thickness of the membrane, the permittivities and conduc- 
tivities of three different regions (i.e. the cytosol, the membrane and the external medium) 
and the properties of the cytosol dispersion. Using Mathematica these parameters can be 
varied individually. Each time only one parameter is varied, while the rest are kept at the 
values in Table I. These variations show interesting behaviours. 

As shown in the figures, there are three sub-dispersions: Aei being the co-field peak 
related to the cytosol, Ae 2 being the anti-field peak related to the membrane and Ae 3 being 
the anti-field peak related to the intrinsic dispersion of the cytosol. 

In Fig. 1, the high-frequency co- field dispersion strength Aei remains relatively constant 
from a e = 1 x 1CT 5 S/m to about a e = 0.01 S/m and then decreases with increasing a e , due 
to a significant reduction in the conductivity contrast between the cytosol and the external 
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medium. Its corresponding characteristic frequency u\ also remains relatively constant in the 
mentioned range and then increases with increasing a e . The anti- field dispersion strengths 
Ae 2 and Ae 3 and their corresponding characteristic frequencies cu 2 and u; 3 show more inter- 
esting behaviours. Ae 2 and Ae 3 swap at between a e = 0.00018 S/m and a e = 0.00019 S/m, 
while u 2 and u 3 show level-repulsion, i.e. their values gain closer, being closest at the same 
value of a e as when the swapping occurs, and then their values move apart again. This 
phenomenon is very common in many physical systems and is frequently observed in atomic 
physics. These interesting phenomena are evidences that both Ae 2 and Ae 3 are real (as op- 
posed to virtual solutions arising from inaccurate calculations) and are common in varying 
many of the parameters, as shown below. 

In Fig. 2, increasing <Xj causes Aei to increase from negative (anti-field) to positive (co- 
field) and then remain constant, Ae 2 to decrease to a constant value and Ae 3 to remain 
constant throughout. lo\ increases monotonically while oj 2 and u; 3 remains roughly constant. 

In Fig. 3, varying uj d has negligible effect on Ae x and thus also Ui. In fact, Ae 2 and Ae 3 
are not very much affected if not for the swapping occurring at about u d = 30000 rad/s. 
Their corresponding characteristic frequencies also show level-repulsion, as in previous cases, 
with the closest point also at about uo d = 30000 rad/s. 

In Fig. 4, Aei and Ae 3 (and also their corresponding characteristic frequencies uj\ and 0U3) 
show negligible variations. Both Ae 2 and u 2 remain relatively constant before increasing. 
They being the only affected ones because they are, as well as the concerned parameter a m , 
related to the membrane. 

Compared with the isotropic mobile charge model with a non-dispersive homogeneous 
core previously investigated, it is interesting that the variations of different permittivities 
and conductivities show remarkably similar results, with the most noticeable difference that 
the swapping and the level-repulsion did not occur in the previous model. 

Fig. 5 shows the real and imaginary parts of the CM factor against the field frequency for 
several values of the medium conductivity, in an attempt to illustrate the results in Fig. 1. 
In this figure, two dispersions are observed. In fact, as shown in Fig. 1, the third dispersion 
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strength is small enough to be neglected, and hence the third dispersion in Fig. 5 cannot be 
shown, as expected. Similarly, we are able to adjust the other parameters respectively, like 
the cytosolic conductivity, circular frequency of cytosol dispersion and external conductivity, 
in order to illustrate the results in Figs. 2~4. However, all of them should show a framework 
similar to Fig. 5, and hence are omitted. 

V. SMALL-GRADIENT EXPANSION 

After investigating models with homogeneous cores, it is natural for us to proceed to 
investigate models with non-homogeneous cores. Here we choose to investigate a model that 
consists of a dispersive core with graded dielectric profile, and a non-dispersive membrane. 
We consider a graded permittivity profile 

e .(r) = €i + ar + 0[a} 2 , < r < R h (30) 

where a is a gradient, which has the unit as permittivity per unit length. We start from 
Eq.(21) and replace the cytosolic permittivity by an equivalent permittivity ej(r) [16] 

€i(r) = e i + ^aR i + 0[a} 2 . (31) 

As shown below, 6j(r) can formally be calculated using the small-gradient expansion of the 
differential effective dipole approximation (DEDA) [16,17]. 

After the substitution, we can expand the CM factor using Taylor's expansion, up to the 
second order: 

U(a) = U(0) +aU'(0) + 0[a} 2 . (32) 

The first term is the same as U iso , while the second term is the correction due to the graded 
profile. 

We now extract the DDSR of the second term as usual, by replacing the permittivities 
by their complex counterparts as in Eqs.(22)-(24), using the following substitution 



10 



X* = — (33) 



Using these substitutions, we can see that the dielectric profile is 



e *( r ) = Q + ^±^ + 0[a] 2 (34) 

1UJ 



Assume the second term has the form 



aU'(O) = °° + ° lW + ° 2W2 + ° 3Vj3 (35) 



C + dw + C 2 w 2 + C 3 w 3 



(36) 



V q + \ €q (17) 

L^i *| I iui 1 /i _|_ icj\2 ^ ' 



((i + u>M)(i + ^M)) 2 

where w = iu and Cs and Ds are constants in terms of the parameters of the model. 

Although the denominator is a quartic equation, there are only two distinct solutions 
for the characteristic frequency. This is due to the differentiation performed in the Taylor's 
expansion, causing each frequency to split into a repeated root. 
By partial fraction aU'(Q) takes the form 

2 Ae 2 A 2 e 

iu> ' ' t-y I iu. 

Wq 9=1 ^ Ulq' 

To avoid confusion, it should be remarked that A 2 e g does not equal the square of Ae q . 
To solve for u q , solve the quartic equation 

(1 + D lW + D 2 w 2 ) 2 = 0. (38) 

u q s are minus the solutions to this equation. 

Aei and A 2 e x in terms of the constants Cs and the characteristic frequencies u q s are 

cji^ 2 (2C - Ci(wi + u 2 ) + wi(C 3 wi(wi - 3^ 2 ) + 2C 2 cj 2 )) 
Aei - , (6V 

[ui - u 2 y 

A 2 £i = (Co - uhjCi + 1M-C2 + C 3 m)) 

1 (W! - CJ 2 ) 2 

Ae 2 and A 2 e 2 can be obtained by replacing oj\ with u 2 and uj 2 with u;i. 

We are now in a position to show how to find ej(r) from DEDA [16,17]. For the dipole 
factor of a graded spherical particle, the following differential equation holds [16,17] 
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% = "3^) K 1 + 2h ^ - ^ - ^^WlKl + + 2(1 - &)*(,•)] (41) 

where b is the dipole factor, r is the radius, and 6j(r) is the dielectric profile. 
Since 

*> = irfe (42) 

solving for b(r) is equivalent to solving for €i{r). 

Since we are doing a small-gradient expansion, b(r) can be expressed as 

b(r) = b + h + 0[a] 2 (43) 

where 

6 » = S (44) 

and 6i can be solved from the differential equation 

dbi 3[(ei + 2e e ) 2 6i - 3e e ar] 



rfr r( Ci + 2e e y 



(45) 



The solution reads 



h ^ M^kr (46) 

using the initial condition at r = 0. 

After putting all the pieces together, e«(r) comes out as in Eq. (31). 

Using the parameters as shown in Table I and a = 0.025/ Ri (where Ri = R e — d and is 
the internal radius), we have done some numerical calculations. This value of a corresponds 
to a change of 10% over the internal radius. The results are shown in Table II. These 
are corrections to the calculations in the previous session due to a small gradient. The 
characteristic frequencies remain the same as in the previous isotropic electrostatic model 
(except now each is a repeated root), while the dispersion strengths are smaller than those 
in the previous model by one to two order of magnitude. This shows that our small-gradient 
expansion is valid. 
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VI. DISCUSSION AND CONCLUSION 



Here a few comments are in order. In view of our recent success in the DDSR of single- 
shell spherical cell model, we are prepared to illustrate the DDSR in various different sit- 
uations. We would like to extend DDSR to cell suspensions of higher concentration. At a 
higher concentration, we expect mutual interactions among cells and the dielectric behaviors 
can change significantly. We may extend DDSR to polydisperse cells, because the cells may 
have different sizes and/or permittivities. The polydispersity can have nontrivial impact on 
their dielectric behaviors. Eventually we have to overcome the analytic continuation, and 
analyze the dispersion spectrum of the full anisotropic mobile charge model and the graded 
cell model. 

Regarding the applicability of the Clausius-Mossotti approach, one can solve the electro- 
static problem first, and then extend to complex permittivities accordingly, as pointed out 
by Jones [1]. As a matter of fact, there are already theories of e.g. Maxwell and Wagner, 
and Rayleigh for heterogeneous dielectrics [18]. In this regard, it is of value to compare these 
theories with the present approach. 

In the present paper, we have discussed isolated particles in the dilute limit. In fact, for 
higher volume fractions, we can use the effective- medium theories instead [10], like Maxwell- 
Garnett approximation or Effective Medium Approximation. 

Throughout the paper, the cells under consideration exist in the form of a spherical 
shape. In this connection, we may include the shape effect as well. More precisely, we 
can extend the graded spheroidal cell model of Huang et al. [17] to include an intrinsic 
dispersion in the core. This is a nontrivial extension and we believe the non-spherical shape 
will have significant impact on the dispersion spectrum. We can consider the following 
agenda: (1) homogeneous spheroidal cell with intrinsic dispersion, without shell [10]; (2) 
graded spheroidal cell with intrinsic dispersion, without shell. Also, item 2 will be studied 
in the small-gradient expansion. 

In view of the present interesting results, the corresponding experiment is suggested to 
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be done. In doing so, one may use coated colloids having a graded core. 

In summary, we have considered a single-shell model with an inhomogeneous graded cy- 
tosol. Realistic cells must be inhomogeneous due to the compartment in the interior of cells. 
In such a model, the cytosol can have a conductivity profile which varies along the radius of 
the cell. A small conductivity-gradient expansion for the DDSR of single-shell graded cell 
model has been done, based on the differential effective dipole approximation [16,17]. We 
have assessed the effects of a conductivity gradient in the cytosol on the dispersion spectrum. 
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TABLES 



Parameters 


Symbols 


Numerical Values 


Cell radius 


Re 


9.5 fim 


Membrane thickness 


d 


8 nm 


External permittivity 


e e 


80e 


External conductivity 


cr e 


1 mS/m 


Cytosolic permittivity 




120e 


Cytosolic conductivity 


<?i 


0.25 S/m 


Cytosolic dielectric increment 


Ae; 


800e 


Membrane permittivity 




7.23e 


Membrane conductivity 


Cm 


4 x l(T 7 S/m 


Circular frequency of cytosol dispersion 


Ud 


10 4 rad/s 



TABLE I. Parameters used for model calculations 



Solutions Symbols Absolute numerical Values 

Characteristic frequencies uj\ 1.03 x 10 8 rad/s 

uj 2 3.18 x 10 4 rad/s 

w 3 1.00 x 10 4 rad/s 



Dielectric dispersion strengths Aei 0.0621 

Ae 2 -0.00133 

Ae 3 1.42 x 10~ 

A 2 ei -0.0613 

A 2 e 2 0.000539 

A 2 e 3 2.02 x 10" 



TABLE II. Results from small- gradient model calculations 
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FIGURES 

FIG. 1. The dispersion strengths (Aei ~ Ae3) and the characteristic frequencies (u>\ ~ ujs) as 
a function of the conductivity of the external medium a e . 

FIG. 2. Same as Fig. 1, but as a function of the conductivity of the cytosol <7j. Typical Oi values 
range from 0.2 to 1 S/m. 

FIG. 3. Same as Fig. 1, but as a function of the circular frequency of the cytosol dispersion uj^. 

FIG. 4. Same as Fig. 1, but as a function of the conductivity of the membrane a m . 

FIG. 5. Real and imaginary parts of the CM factor as a function of the circular frequency of 
the external field. Re[- • •] (Im[- • •]) denotes the real (imaginary) part of • • • . 
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The effect of varying o. on Ae 
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The effect of varying co on Ae 
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